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A complete study for the production of neutral (ft , H°, A°(— 0°)) and charged Higgs (H ± ) 
bosons at electron-photon colliders is presented in the context of the minimal supersymmetric stan- 
dard model. A particular choice of the non-linear _R^-gauge is used to evaluate the amplitudes of the 
reaction ej — > e<^'\ The resulting cross section indicates that it will be possible to detect a signal 
from the neutral Higgs bosons for most regions of parameter space at the future linear colliders with 
' -y/s = 500 GeV through the reaction e7 — > ecf)® . This reaction also offers the interesting possibility 

to measure the Higgs mass through the detection of the outgoing electron. The production of the 
charged Higgs boson (H + ) through the reaction ej — * f e H ± has in general smaller values for the 
cross section, which seems more difficult to observe. 
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' I. INTRODUCTION 

> 

The search for Higgs bosons at future colliders has become the focus of extensive studies, because of its importance 
as a test of the mechanism of electroweak symmetry breaking |l| . The detection of the full spectrum of scalars seems 
Cn| ' necessary in order to determine the nature of the physics that lies beyond the standard model (SM). Among the 
extensions of the SM, supersymmetry has received increasing attention, not only because of its aesthetical properties 
as a field theory, but also because the naturalness problem of the SM can be alleviated within the so-called minimal 
supersymmetric standard model (MSSM) j2| . 

The MSSM contains two Higgs doublets, whose physical spectrum includes a charged pair (H^), two neutral 
Qh! CP-even scalars (h°, H°), and one pseudoscalar (A ). The MSSM Higgs sector is determined at tree-level by two 
parameters, which are nowadays chosen as the A mass and tan/3 (the ratio of the vevs. of the two Higgs doublets). 
| This in turn fixes the values of the neutral Higgs mixing angle a and the remaining Higgs masses, which obey the 
tree-level relations, m^o < mz < m^o , m^o < m^o < m^o , m\y < m H ± . However, these relations are substantially 
modified when the effect of radiative corrections is included |] ||. In particular, it makes possible that mz < rn^o, 
m^o < m h o, and for some regions of parameters m^o could even reach a value of about 130 GeV |^|^]. 

Current LEP2 limits on the Higgs masses are of about 90 GeV for the light scalar (h°), and, depending on the 
value of tan/3, up to about 85 GeV for A j9|. However LEP2 will be able to cover the region up to m^o < 110 
GeV |fio| . Recently, it was found that Tevatron can be used to test a significant portion of parameter space through 
the reaction pp -> Wh° + X |ll[, moreover the possibility to perform b-tagging with a high efficiency has opened the 
window to detect the mode pp — > bbh a + X at LHC |l^] and it will allow to test the large tan (3 region of parameters 
at Tevatron jL3|. On the other hand, it has been shown Jl4| that LEP2 plus LHC results will be able to cover almost 
all the MSSM Higgs sector parameter space, and through a combination of the reactions pp — > it + h°(-^ 77) + X, 
pp — * h°(— » ZZ*) + X |15-|l7|l and weak boson fusion |Ti|, the full region will be covered. At future linear colliders 
like NLC TESLA [2C] or JLC it will be possible to search for neutral and charged Higgs bosons, through 
the production reactions e+e" h°Z, h°^f HIH and e+e ~ ~* H+H- |||. 

The future linear colliders can also operate in the e7 mod e [p5| ; in this case the production of the SM Higgs 
boson has been studied through the 77 — * (jp mechanism prj| , [27| , and also with the full 1-loop two-body reaction 
e7 — > |p8| p0| . The production of the pseudoscalar Higgs boson A through the reaction e7 — > A°e has been 
studied too |[31|| , using the photon pole approximation. Even though, these reactions can occur at tree-level, they 
receive the main contributions at one-loop level where the heavy particles of the model play an important role. These 
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processes could be used to measure the couplings: ^7*7, 0°Z*7 |2jj|^,|33|] , which constitute important one-loop 
predictions of the theory, and are also sensitive to the effects of new physics. Moreover, because of the possibility to 
measure the electron momenta, this reaction offers the possibility to determine the Higgs mass with high precision. 

In this paper we present a detailed analysis of the production of the neutral and charged Higgs of the MSSM at 
7e~ colliders through the reactions e~f — ► e§\ (<$ = h° , H° , A ) and — > vH^, at the one-loop level. Our goal is to 
determine the regions of parameter space where a signal is detectable, and also to find out where it will be possible to 
distinguish between the MSSM and SM Higgs signals. We shall assume that the super-partners are heavy, and thus 
decouple from the amplitudes J34| , however the effect of squarks will be included in the Higgs effective potential using 
the approximations presented in Ref. (3^] . 

The organization of the paper is as follows: Section 2 discusses the production of neutral Higgs bosons. It includes 
a discussion of the region of parameter space that can be excluded, and also on the determination of the Higgs mass. 
Section 3 is devoted to the production of the charged Higgs, whereas the conclusions of our work are presented in 
Section 4. Details about the non-linear gauge used in the evaluation of the 1-loop amplitudes as well as the explicit 
formulae for the various amplitudes are presented in the Appendices. 

II. THE NEUTRAL HIGGSES 

We now proceed to present the results of our calculation of the amplitudes of the reactions: 

1 (k 1 ) + e-( Pl )^cj ) (k 2 ) + e-(p 2 ), (1) 

where <jp denotes any of h°, H°, A , and we have also displayed the notation for incoming and outgoing momenta. We 
have organized the calculation according to the U eBl (l) gauge invariance, thus the diagrams are grouped as follows: 

1. three-point diagrams characterized by the </>°7*7 and by the (j>°Z*j coupling (Fig. 1-a) ; 

2. Z and W-mediate box diagram and its related triangle graphs (Fig. l-b,c,d), 

The triangle graphs related to the Z- and VF-mcdiated box diagrams are the one-loop (f>°e*e 3-point functions. 
These groups of diagrams are finite and gauge invariant by themselves. The group (1) receives contributions from 
loops of charged fermions, W gauge boson, and charged Higgs boson H ± . On the other hand, the group (2) is sensitive 
to the (jPW + W~ and <f>°ZZ vertices, respectively. For A there are contributions from fermionic triangles only, thus 
only group (1) appear. 

In the linear gauge, the reaction — > e<^ receives also contributions coming from the reducible diagrams with 
the Z*~/ self-energy, which is an extra complication because one needs to perform a renormalization of this term. 
However, in the nonlinear gauge this term is absent as consequence of the explicit U cm (l) gauge symmetry in the 
W, goldstone and ghost sectors. Notice that there are no contributions coming from the charged Higgs boson to the 
group (2). This happens because we use the approximation m e = and also because there is no ^yW^H^ vertex at 
tree-level .[] 

A. Production of h°,H° 

The Mandelstam variables used in this calculation are defined by s = (hi +P1) 2 , t = (Jfci — fc 2 ) 2 , and u — [k\ — P2) 2 ■ 
In addition, e A '(fc 1 ,Ai) will denote the 7 polarization vector. We have evaluated the amplitudes using dimensional 
regularization, with the help of the program FeynCalc |3(| and the numerical libraries FF p7| , p8[ . 

The result for the total amplitude of the reaction — > eh°(H°) can be written as: 

M =M y + Mz+M h z OK + M^, (2) 

where M^, Mz, M^ ox , and M^ x correspond to the sets of diagrams (1), (2), (3), and (4), respectively. The 
amplitudes coming from the </>°7*7 and <p°Z*j couplings can be written as follows: 



1 We could think on possible contributions arising from reducible diagrams that include the 7*0° and Z*<p° self-energies (with 
the virtual fields tied to the electronic line). However, it can be shown that these terms always vanish. 
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M i> z = 434 u(p 2 )-f v (a lyZ - fo 7 ,z75)w(pi)e A '(fci, Ai).F 7l z(fci ■ k 2 g^ v - k 2f ,k lL ,), 



(3) 



where 



(4) 
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(5) 



with a 7 = 1, 6 7 = 0, a z = 1 — 4s 2 ,, 6 Z = 1, t„ 



-, and c 2?I 



sj,. The coefficients v^, and characterize 



the 0°//, and cjP^H* couplings in the MSSM, and are given by: 



4^, H°uu 

sin p 5 _ 

COS p 1 

sin /3 ' 

cos p 1 1 



cos(/3- a), 0° = H° 
v<i> ~ i sin(/?-a), 0° = ft , 



cos(/3-a) - 2 i s-cos2/3cos(/3 + a), c/) = H° 
sin(/3 — a) + ^3- cos 2/3 sin(/3 + a) , 



b° = h°. 



The functions f)' 2 , F°, and F\ arising from fermion, scalar and gauge boson loops are given in the Appendix B. 

The amplitude for the contributions of the Z-mediated box diagram, including the related (f) e*~e~ triangle graphs, 
is the following: 



M h z ox = l -^-p^u(p 2 )Y (a z 75 ) 2 w(pi)e^(fci, Ai) 



x [-A(t, s,u) (fci ■ pig^u - Pifik lu ) + A(t,u,s)(ki ■ p 2 g^ ~ p 2 ^k lv )] , 



(6) 



where the functions A(t, s, u) and A(t, u, s) are presented in the Appendix B. 

The amplitude for the VF-mcdiatcd box diagrams and its related triangle graphs is given by: 



M h ^^ l -^^u(p 2 )^(l- l5 ) 2 u(PiV(ki,\i){A 12 (t,s,u) 
z ±s w c w 

x (fci • PWn V - P\pk\v) - A 2 i(t,s,u)(k- L -p 2 g^-p 2f ik lu )} , 



(7) 



where the functions A4 and A^ are given in the Appendix B. 

After squaring the total amplitude, the corresponding cross-section for this process is given by: 



1 
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(8) 



where 
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w 



(-t) [(s 2 + u 2 ) F lZ + s 2 F s + u 2 F u ] , 



(9) 



with 
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\F 1 \ 2 +a\F z \ 2 + 2a z Re(F z F;), (10) 



F s = \A X2 (t, s, u)\ 2 + c\A(t, s, u)\ 2 + 2Re {[A 12 (t, s, u) - 2aA(t, s, u)]F* 

+ [4bA 12 (t, s, u) - 2dA (t, s, u)] F* - 2A 12 (t, s, u) A* (t, s, u)} , (11) 

F u = F s (s^u), (12) 

where we have defined a = 1 + a 2 , b = 1 + a z , c = 1 + 6a 2 . + a*, and d — a z (a 2 + 3). 

Finally, in order to obtain the total cross-section (ot), one needs to convolute a with the photon distribution, 
namely: 

~ ~q I f^HWs, (13) 



where 5 denotes the squared cm. energy of the e + e~-system, and the photon distribution is given by: 

1- x + 



/-,= 1 



D(0 
where 



1 4x 4a; 2 

(14) 



1-x £(l-x) £ 2 {l-x) 2 



w-^-rfjw+Q+l+rwhr- (15) 

Notice that, as in Ref. [p9| , one should use, instead of the photon distribution, an exactly monochromatic initial 
photon beam which can help in distinguish the phisical effects related to the particular collision process from details 
depending on the final realization of the laser beam. 

To evaluate the cross-section for h° and H°, we have taken m f = 175 GeV, and the values for the electroweak 
parameters given in the table of particle properties (3^]. In fact, as it was discussed in Ref. p3], the contributions of 
the boxes can be neglected. The cross section for h° is about 4.2 fb for m^o > 200 GeV and E c . m . = 500 GeV, as can 
be seen from Fig. 2. With the expected luminosity at the future linear colliders [fl9|-f2l| of about 50 fb~ 1 /yr, it will 
be possible to observe up to about 210 h° + e events, which should allow to study the properties of the Higgs boson. 
Among them, it will be interesting to test its spin and CP-even nature by studying the angular distributions, however 
this aspect is beyond the goal of the present work, where we are mainly interested in determining if the MSSM Higgs 
sector can be tested at electron-photon colliders. Total cross-sections for the production of the heavy Higgs boson 
H° as a function of the pseudoscalar mass m^o and for different values of tan/3 are shown in Fig. 3. These graphs 
present the total cross-section as a function of the pseudoscalar mass (m^o) for several values of tan p for cm. energy 
of 500 GeV. 



On the other hand the amplitude for the production of the pseudoscalar contains only contributions from the 
fcrmionic triangles. Thus, the amplitude takes the form: 



where 



M=M 1 +M Z , (16) 

-ia 2 Q 2 N c fpm w . 4m? , 2 2 s ( c^pkfk^ \ 

M 7 = L u(p 2 )-f y M(pi)e' i (fci, Ai)— ^-Cq (t,m%mj) ' 1 , (17) 



and 



ia 2 Q } C } v N c fi3m w Am) ( e^ a (}kfk^ \ 
Mz = —. 3 2 u{p 2 )j (C£ - C^7 5 )tt(pi)e' 1 (fei,Ai)— j-C [t, m% mj) — , (18) 
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with Cn ( t, m 



'o(t 



C, 



o(V. 



i A , 0, rrif, m 2 , m 2 



the Passarino-Veltman three-point scalar function written in the 



notation of the FeynCalc program and 



cot /?, / = u 
tan /3, f = d, 



The total squared amplitude is 



-t s 2 + u 2 
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f 2 °v 2 
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f 2 



t 2 



s~ie (if 



2 s 2 r 2 



Q f t- 



16si 



c wQ 2 f 



\2m 2 f Co (t, m 2 A , m 2 ) 



(19) 

(20) 
(21) 



To obtain the total cross-section we use the expressions written previously for the CP-even Higgs bosons (Eqs. 13-15). 
Results are shown in Fig. 4, where we plot again the total cross-section as a function of the pseudoscalar mass (m A a ) , 
for several values of tan/3 and for cm. energy of 500 GeV. It can be noted that the cross-section for A a is smaller 
than the one resulting for h°/H°. Thus, in this case it will be more difficult to detect the signal. 



C. Backgrounds and exclusion contours 

The final signature for the reaction ej — ► ec/)® depends on the decay of the Higgs boson. For h , the dominant mode 
is into bb, whereas for H° and ^4° this decay can also be relevant for some regions of parameter space. To evaluate the 
B.R. of the Higgs bosons we shall assume that the decays into SUSY modes (i.e. charginos, neutralinos, sfermions) 
are not allowed and take the relevant equation for the decay widths from H. Our results are in agreement with the 
ones obtained in the literature, for instance in ref. fl40| . Thus, we shall concentrate on the signature coming from the 
decay 0? — * bb. In this case the main background comes from ej — ► ebb, which receives contribution from 8 graphs 
at tree level; we have evaluated numerically this processes using CompHep ^l],|^]. Following Dicus et al. |pl|] , we 
have also imposed a cut on the angular distribution of the outgoing electron (| cos#| < 0.98, relative to the incident 
photon) that reduces significantly the background, while retaining most of the signal rate. To determine the region 
of parameter space, which is taken as the plane tan/3-m^o, where the Higgs signal is detectable, we have proceeded 
as follows: for each value of m A a we evaluate the Higgs masses and the cross-sections, then we find the value for tan /3 
where the cross-section of the signal is above the background at the 3 and 5a level. 

In Fig. 5, we show the regions in the plane tan/J-m^o where the cross-section coming from h° has detectable 
values. The region to the right from the heavy line is where the signal is detectable, above backgrounds at the 3tr 
level. The dashed line denotes the contour at the 5tr level. It can be seen that the signal is above the backgrounds 
for a significant region of parameter space. 

On the other hand, Fig. 6 shows the region where the cross-section from H° reaches detectable values. Again, 
The region to the left from the heavy (dashed) line is where the signal is detectable at the 3 (5) a level. It can be 
noticed that this region covers the sector of parameter space where it is more difficult to detect the light h°. Thus, the 
cross-sections for h°/A° and H° /A play a complementary role in providing a detectable signal for the full parameter 
space. Moreover, we also notice from the superposition of Figs. 5 and 6, that there is a small region where the two 
Higgs boson signals can be detected, which will allow to distinguish clearly between the MSSM and the SM Higgs 
sectors. 

Finally, we also want to stress the fact that this reaction offers a unique opportunity to obtain a clean measurement 
of the Higgs mass, thanks to the possibility to measure the outgoing electron momentum. The mass of the Higgs 
boson is related to the maximal energy of the electron as: 

m h a =Vs— E max . (22) 

Thus, the precision attained for the electron energy will translate into a good determination of the Higgs mass. 

On the other hand, the largest values of the cross-section for ^4° are obtained for large values of tan (3 (which is 
usually assumed to be at most of order 50); however, we find that even in this case the resulting cross-section does 
not seem to give a detectable signal. Moreover, for large values of tan/3 there appears a mass-degeneracy between 
A and h° or H°, which will make difficult to distinguish the individual signals. In this case, it will be neccesary to 
optimize the cuts to be able to detect the pseudoscalar A , and to separate it from the largest signals coming from h° 
and H°. Otherwise, one would have to add the respective signals; for instance, one could imposed the criteria that 
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whenever the difference between the Higgs masses is less than 5 GeV (a conservative criteria given the possibility to 
obtain a precise measurement of the Higgs mass), the individual contributions to the signal must be combined. We 
find that this only helps to enlarge the exclusion contours for tan {3 above 30. 

Another option to detect A (and a heavy H° too) is to use the planned second and third stages of the future 
e-gamma colliders -^T[ , which could reach an energy of 1 TeV and 2 TeV, with integrated luminosities that could 
reach 125 and 500 fb~, respectively. 



III. THE CHARGED HIGGS 



Now, turning to the production of the charged Higgs, we observe that it can also proceed through the 1-loop 
reaction: 



7(fci) + e (pi)-»H (*a) + f e (pa). 



(23) 



The diagrams encountered in the calculation of the e~7 — > vH~ are shown in Fig. ||. They include: triangle graphs 
with bosons and fermions in the loop (Fig- §k), seagull-type graphs with bosonic contributions (Fig. ||b-c), and finally 
those with self-energy insertions (Fig. gd, [83, |8|f). Fig. ^f gives a vanishing contribution for massless external fermions. 



A. Production of H 



Our result for the total amplitude is written as: 



M 



^2)7^(1- 75 MPiK^Ai) 



x (Vf + V H+r + V W4 ,o) (k 2ll k 1 k x ■ k 2 ) + iAfe^apk^kt, 



(24) 



where Vf, Vjj+^o, Vw^o and Af denote the contribution from the different sets of graphs shown in Fig. |S| and are 
given in the Appendix C. 

Fig. [9] shows the cross-section for our process for several values of tan/3 (= 1.5, 5, 10). The cross section is small, 
about 0.5 fb for m H ± = 200 GeV and E c m . = 500 GeV, which can give 25 events with an expected future linear 
colliders luminosities of 50 fb _1 /yr. In Fig. |To| we compare the cross section from e-f — > H~i>, with the pair production 
e + e~ — * H + H~, c~f — > H + H~e and 77 — > H + H~ . The production mechanisms e + e~ — * H + H~, 77 — * H + H~, 
have been discussed in the literature [j43[ . On the other hand, the reaction — ► H + H~e is evaluated using the 
Williams- Weizsacker approximation (for the second photon) and we use the sub-reaction 77 — > H + H~ . It can be 
seen that the single production dominates only for large Higgs masses i.e. for values that lay beyond the threshold 
for pair production. 



B. Backgrounds 

In this paper we only show (Fig. ^) the results for fixed values of tan/3 = 1.5, 5, 10, and as can be appreciated it 
turns out that the resulting cross-section has smaller results, which hardly seem detectable. In order to determine 
if the signal is detectable, one needs to consider the potential backgrounds, which depend on the Higgs mass, since 
this determines the decay signatures. For instance, if m#± < m t + m&, the dominant decay mode is H + — * tv 
which reaches a branching ratio of order 1 for most values of tan/3. In this case the background will come from 
the production 7e — > W*v, which for masses m H ± « mw will be much larger than the signal (it reaches a ~ 4pb) 
and probably will not allow detection. For masses somehow larger, both the signal and the background will be more 
suppressed and the question of detectability will depend on the experimental ability to identify the decay mode and 
to reconstruct the charged Higgs mass. 

For heavier Higgs masses, M H ± > m t + m&, the dominant Higgs decay is into t + b, and in this case one needs to 
compare the signal with the background arising from single top production, e"/ — * vtb for M H ± > m t + mb, which 
according to the results of Boos et al. p4| , has a ~ 15fb for ^/s = 500 GeV and m t — 175GeV. Thus we can see that 
at the level of total cross-section the background is again larger than the signal. However if one makes a cut in the 
invariant mass of the t — b system, then it will be possible to reduce the background. 
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IV. CONCLUSIONS 



We have studied the production of the neutral CP-even and charged Higgs boson of the MSSM (h° , H° , H^) at 
future e7 colliders, through the reactions ej — ► eh , eH°, eA°, uH~ . The amplitudes are evaluated using a non-linear 
i?^-gauge, which greatly simplifies the calculation. The resulting cross section indicates that it is possible to detect 
the light neutral Higgs boson (h°) for most values of parameters. On the other hand, detection of the heavy neutral 
Higgs bosons H° seems possible only for light values of m^o . We have determined the regions in the plane tan 
where the signal is above the backgrounds. The cross-sections for h° and H° play a complementary role, since the 
region where H° reaches detectable values occurs precisely in the region where it is more difficult to detect the light 
h°. Thus, both reactions allow to cover the full plane tan/3 — m^o with at least one detectable signal. However, it is 
found that the possibility to distinguish the MSSM from the SM case, through detection of both h° and H° signals, 
occurs only for a limited region of parameter space. 

On the other hand, the results for the pseudoscalar and the charged Higgs boson have smaller values of the cross- 
section, which seems difficult to detect. 
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APPENDIX A: THE NON-LINEAR GAUGE 

In the evaluation of the above processes, we have used a Feynman-t'Hooft version of the nonlinear gauge |45|-|49|l 
that greatly simplify the calculation. In order to appreciate the advantages that offer a nonlinear gauge, we find 
convenient to compare it with the conventional linear gauge for the SM. It is well known that in a linear gauge 
only the charged fermion sector presents explicit electromagnetic gauge symmetry (U em (l)), since the gauge fixing 
procedure used for the SU(2) sector destroys manifest U em (l) symmetry in the charged gauge boson (W ± ) and ghosts 
sectors (c, c), when they are considered separately, because the gauge functional used to define the H^-propagator does 
not transform covariantly under the U em (l) symmetry. It follows that instead of obeying naive Ward identities, these 
sectors are related through Ward-Slavnov identities. In order to obtain a finite and U om (l)-invariant result for higher- 
order (loop) calculations, one must sum over the contributions arising from the W gauge boson, the charged ghosts, 
and the W A± G T combined effects, denotes the would-be Goldstone boson. On the other hand, the functional 
used to define the PF-propagator in the non-linear gauge contains the electromagnetic covariant derivative. Thus, the 
Ucm(l) symmetry is respected by each charged sector of the SM. It follows that a finite and U em (l) gauge invariant 
result is obtained for each type of diagram containing a given kind of charged particles. Moreover, the number of 
diagrams involved is considerably reduced because there are no M /± G =F combined effects. In the SM, the functionals 
that define the W, Z, and 7 propagators are given by: 

/+ = - itm w G + , 

f z = d^Z" - tmzG", (Al) 
f A = M". 

where — — ig' B^, with B^ = —s w Z^ + CwA^, G° is the neutral would-be Goldstone. £ is the gauge parameter, 
which in general is different for each of the three gauge bosons. However, in the non-linear Feynman-t'Hooft version 
of this gauge, one takes £ = 1 for all sectors. The gauge fixing Lagrangian is given by: 

c G F = -\f + r-^{f z ?-^{f A )\ (A2) 

which in turn removes the H /± G =F 7 and W^G^Z vertices from the Higgs kinetic energy term. Notice that the f + f~ 
term is U em (l) gauge invariant, which implies that the Fadeev-Popov Cpp Lagrangian is also U cm (l)-invariant. It is 
also possible to proceed further with this scheme and remove more un-physical vertices, such as <jPW G^^f{Z) and 
G°W ± GT 7 (Z). This procedure requires to define the functionals nonlinearly both in the vector and scalar parts. 
However, for the present purposes it is sufficient to use the scheme presented above. 

In this case only the couplings in the gauge and ghost sectors are modified. We shall give the lagrangian for these 
sectors in a form that clarifies the role of the electromagnetic gauge invariance. For this, we define the derivative: 
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Dp = — igW^, with = c w Z^ + s w A^. This operator shares with the derivative D defined above the property 
of containing the cm. covariant derivative. After adding the gauge-fixing lagrangian, one obtains the couplings for 
the gauge bosons (W, Z, A), which are contained in: 

C = -\{D^W+ - b u W+)\D»W +lJ - D»W+n - ig(s w F^ + ^Z^W^W^ 



2 

-ig 



s w F„ v + c w Z^ + i 9 - {W-W+ - W~W+) 



w-»w +v 



-\Z, V Z^ + \m\Z.Z* 1 (d,Zn 2 -F„F, V - 1 (d^f (A3) 

where F^ u — d^Av — d„A^, Z^ v = d^Zv — d^Z^. Notice that, apart from the gauge fixing term d^A^, this lagrangian 
is U em (l) invariant. 

On the other hand, the ghost sector is substantially modified in the non-linear gauge; however, since the gauge- fixing 
functionals are U(l) invariant, the charged part of this sector is also invariant. We present now the corresponding 
Fadeev-Popov lagrangian, written in such a way that U om (l) invariance is explicit, namely: 

Cfpg = -c- [d^D" + £m w (m w + 0° + zG°) 

2 

Cw 

-igc w W+^ (d^c z ) c- - ieW + » (d^) c~ 

—gc2w(,mzG + (c~cz + czc~) — e£mwG + (c"c 7 + c 7 c _ ) + h.c. 

-c z [□ + im z (m z + 4> )] c z - c 7 Dc 7 (A4) 

where c ± (c ± ),cz(cz),c 1 (c 1 ) denote the pairs of ghosts associated with the W 7 Z,A gauge bosons respectively. The 
phase-convention for the charged ghost is (c+)t = c~ and (c) 1 * = cT . In addition, one has: (fP = cos(a — /3)H° + 
sm(a — where h°,H° denote the light and heavy neutral Higgs bosons of the MSSM. Notice that neither the 

charged Higgs (H^) nor the CP-odd Higgs bosons (A ) appear in this lagrangian, which can be understood if one 
reminds that these fields do not appear in the definition of the gauge fixing terms, and also because they do not 
receive a v.e.v.. 

Finally, we would like to mention the properties of the couplings of gauge bosons with the higgs and goldstonc 
bosons. First, in this gauge the couplings W + G~'f and W + G~Z are absent by construction, whereas the couplings 
G + G~ Z(j), H + H~ Z(j), <f) a W + W~ and W + H~(fP do not depend on the gauge-fixing procedure. However, the cou- 
plings (f)°G + G~ or (f)°c + c~ can depend on the choice of the gauge-fixing terms, when both the vector and scalars 
sectors are chosen non-linearly. 

In conclusion, the Feynman rules needed for our calculation, which are depending of the gauge fixing procedure, 
arise from the lagrangians of Eqs. (A3-A4) written above. 



APPENDIX B: LOOP FUNCTIONS ASSOCIATED TO NEUTRAL SCALARS 

In this appendix we present the various loop functions coming from the e~7 — ► e"(fp processes. The functions 
arising from fermion, scalar and gauge boson loops, which characterize the </>°7*7 and (f> a Z*j couplings, are given by: 

4 777 C 1 

if 2 = 1 1 - \ {ml - t - 4m}) Co (t, m\, mj) 

+ ^T— t [Bo K, rn}) B a (t, m))} j , (Bl) 

4m 2 i t I 

F ° = -Z^7 1 + 2 <± C o {t,m%m% ± ) + -3— [B (m%m 2 H± ) - B (t,m 2 H± )] , (B2) 
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F 1 

7 



4m^ 



2771^ 



3 1 



2771^ 



2t 



"w 



2m^C (i,m|,m^) 



+ 3 + 



2m^V 



mi - t 



(B3) 



Fi = -4 (3 - 4) m 2 w C (t, m% m 2 w ) + -Hj^L 



5 + 



2m^ 



1 + 



2771^ 



x I 1 + 2m^C (t, m\, m^) H 3 1 [ B o (m|, m^,) - B (t, m^)] } , 



,2 „,2\ 



(B4) 



where C (i,m|,m?) = Co(0,i,ra^,raf,m?,ra?), S (m^,m?) = B (m^,m?,mf) and B {t,m 2 ) = B (t,m 2 ,m 2 ) are 
Passarino-Veltman three- and two-point scalar functions written in the notation of the FeynCalc program. 

The loop functions coming from Z-mediated box diagram, including the related </>°e*~e~ triangle graphs, are given 
by: 

4(t, a, u) = ^ | [m|(s + «) - «u] D (l, 2, 3, 4) + (a - m|) [C (l, 2, 4) 



-C (l, 2, 3) - — C (2, 3, 4) + - ( t - s 
s s s 



2t 



2m 2 z st 



(s + t){s-m 2 z ) 



C (l,3,4) 



+ -— [B (3,4)-B (l,3)] 
s + i 



(B5) 



u, s) is obtained from A(t, s, u) by means of the interchange s <-> u. The arguments of the scalar functions arc 
given by: 

D (l, 2, 3, 4) = D Q (0, s, m|, u, 0, 0, m^, m^, m|, m|) , 
C (l,2,4) = Co(0,s,0,m2,m2,m|), 
C (l, 2, 3) = Co (0, 0, u, m*, m 2 ? m |) ; 
C (2,3,4) = C (s,0,ml,m%,ml,m%) , 
C (l,3,4) = C" (0,m|,«,m^,m|,m|) , 

73 (3,4) = B (m|,m|,m|), 

B (1,3)=B (u,m*,m|). 



(B6) 



Finally, the loop functions associated to the T'F-mcdiated box diagrams and its related triangle graphs are the 
following: 



A 12 {t, s, u) = 2<4 [A^t, s, it) + A 2 (t, u, s)] , 
A 21 (t, s,u) = 2c* [A 2 (i, s,u) + i4i(t,u,s)] , 



(B7) 



with 



Ai(t,s,u) = 



1 [ s — m 



2st 



+ (s - m^) 



,r [m^(s + u) + st] D (l,2,3,4) 



C (2, 3, 4) - -C (l, 3, 4) - IZJf C (l, 2, 3) 



+ !i±*C (l,2,4) 
s 



2f 

+ — — [B (l,2)-B (l,3)] 
s + u 



(B8) 



A 2 (i, s,it) = 



t + u 



m 



w 



2tu 



(m^(s + u) - st) - 2rr? w t 



73 (1,2,3,4) 
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+ (t + u - 



in 



w 



-C (2,3,4) + -C (l,3,4) 
u u 



u 
2t 



-C (l,2,3) + 



u 



2ut - t" 



u(t + u) 



Co(l,2,4) 



t + u 



[Bb(2,4)-B (l,2)] 



2t 



s + u 



[So(l,3)-B (l,2)] 



} 



(B9) 



The scalar functions have the following arguments: 



A>(1,2,3,4) 
Co(2,3,4) 
C (l,3,4) 
0,(1,2,3) 
Co(l,2,4) 
So(l,2) 
Bo(l,3) 
So(2,4) 



Do (0, 0, 0, m|, i, s, m^, 0, m^, m^) , 
Co (0,0, s,0 

Co (0,0,t,m^,0,m^) , 



C o (t,0, 



m^m 



0. 



Co (0,s,TO^,m^,0,m^) , 
S (m|,m^,m^) , 
-Bo (i, m w) . 

B o (*,0,m^). (BIO) 

The expressions for Ai(£, u, s) and A 2 (i,u, s) can be obtained from the respective Ai(£, s,u) and A 2 (i, s,w) through 
the interchange s <-> m. 



APPENDIX C: LOOP FUNCTIONS ASSOCIATED TO CHARGED SCALARS 



In this appendix we shall write the explicit formulae for the functions Vf, Vjyi^o, Vw<jfi and A/ using a notation 
that, hopefully, facilitates its use by the interested reader. Our convention for the scalar functions Bo, Co is the 
same as in the previous Appendix. In this case, the contributions arising from loops with fermions, scalars and gauge 
bosons, have the following expression: 

V f = J2 { tQu + Qd) (m 2 d tan 8 + m 2 u cot 0) 

f 1 m H± 

+ 2Qd [ml (ml tan (3 + m 2 u cot 0) - (t - m 2 H± ) m 2 d tan 8] C (0, m 2 H± ,t,m 2 d ,m 2 d , m 2 u ) 
+ 2Q U [ml (ml tan/3 + m 2 u cot 8) — (t — m 2 H ±) m 2 u cot 8] C (0, m 2 H ± , t, m 2 u , m 2 u , ml) 

t2 i„„ a i ^,2 



and 



+ 



m 2 , tan 8 + mf. cot , , , , _ , 9 N 
d , 2 " (m 2 u -m 2 d + (Q u + Q d ) m 2 H± ) 

~ m H ± 



+ 2 (Q dm d tBxi (3 + Q u m u cot 0)\ [B (t, ml, m u ) - B (m H ± , m d , mj] 



+ m u md- 



m 



[B (m|± , m^, m£) - B (0, m^, m£)] 



4f = X] ^ {Qd m d tan /3C (0,m 2 H± ,t,m 2 d ,m 2 d ,ml) - Q„ra^ cot /3C (0,m 2 H± ,t,m 2 u ,m 2 u ,m 2 d ) 
f 

ml tan /3 + cot 8 



+ 



t - m 2 H± 



[B (t, m 2 d ,m 2 u ) -B (m 2 H± ,m 2 d ,ml 



(CI) 



(C2) 



( m 2 - 



t - m 2 H± 



{l + 2m^±C (0,t,m 2 H± ,m 2 H± ,m 2 H± ,mlo) 



+ 



+ 



m 2 H± - 1 



t-m 2 H± 



[B (t, m 2 H± , m\ ) - Bo (m 2 H ± , m 2 H± , m| )] 



[B (0, m^± , m^ ) - B (m^± , m 2 H± , m| )] > , 



(C3) 
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m lo + mi 




+ (3m 2 H± - m\ + m 2 w - 2t) C (0, m 2 H ± , t, m^, m^r, m^o) 
(m 2 H± - m 2 - m 2 ^ (m 2 o - m 2 w - tj - 4m^ (m 2 H± - 1 



% m w {p ~ "' 



2 



V " / V H±~ m l 0+m w) „ „ 2 2 2 M I 

H ^2 2 l B v (°> %> TO 0° ) _ B o {m H ± , /%, to 0O J J > , (C4) 

zm w m H± 



(m 2 - m^) ( 
2m 

where cr^o = —Vh<>, a h a — —v^o and 4>\ = H , 0° = ft 
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FIG. 1. Classification of graphs that contributes to the reaction cy — > eh°. 



FIG. 2. Total cross section for the reaction ej — > eh° for different values of tan /3, with yfs = 500 GeV. 



FIG. 3. Total cross section for the reaction ej —> eH° for different values of tan/3, with yfs = 500 GeV. 



FIG. 4. Total cross section for the reaction cy — > eA° for different values of tan /3, with yfs — 500 GeV GeV. 
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FIG. 5. Cross section contour lines for the reaction ej — > eh°, for different values of ff(fb), with y/s = 500 GeV. The region 
to the right from the heavy (dashed) line is where the signal is detectable (above backgrounds at the 3 (5) a level). 



FIG. 6. Cross section contour lines for the reaction ej — > eH°, for different values of cr(fb), with yfs = 500 GeV. The region 
to the left from the heavy (dashed) line is where the signal is detectable (above backgrounds at the 3 (5)<r level). 



FIG. 7. Cross section contour lines for the reaction ej — > eA°, for different values of cr(fb), with ^fs = 500 GeV. 



FIG. 8. Classification of graphs that contributes to the reaction e 7 — > v e H . 



FIG. 9. Cross section contour lines for the reaction e 7 — > vH , for different values of tan (3, with y^s = 500 GeV. 



FIG. 10. Cross section for the reaction e 7 — > vH , with y/s — 500 GeV compared with different production channels. 
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